Application of Gradient-Corrected Density 
Functional Theory to the Structures and 
Thermochemistries of SCF3, TiF4, VF5, and 



Density functional theory(DFT) and Hartree-Fock(HF) calcula- 
tions are reported for the family of transition metal fluorides SCF3, 
TiF4, VF5, and CrFg. Both HF and the local-density- aproximation 
(LDA) yield excellent agreement with experimental bond lengths, while 
the B-LYP gradient-corrected density functional gives bond lengths 
0.04 — 0.05 A too long. An investigation of various combinations of 
exchange and correlation functionals shows that, for this series, the 
origin of this behavior lies in the Becke exchange functional. Much 
improved bond distances are found using the hybrid HF/DFT func- 
tional advocated by Becke. This approximation also leads to much 
improved thermochemistries. The LDA overestimates average bond 
energies in this series by 30 — 40 kcal/mol, whereas the B-LYP func- 
tional overbinds by only ~ 8 — 12 kcal/mol, and the hybrid HF/DFT 
method overbinds by only ~ 2 kcal/mol. The hybrid method predicts 
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Abstract 



the octahedral isomer of CtFq to be more stable than the trigonal 
prismatic form by 14 kcal/mol. Comparison of theoretical vibrational 
frequencies with experiment supports the assignment of an octahedral 
geometry. 
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1 Introduction 



The new generation of "gradient corrected" density functional have been 
shown to yield remarkably accurate thermochemistries for organic systems. 1,2 
Applications to inorganic complexes are less extensive, but those studies 
which have been reported are very encouraging. 3-7 We have recently devel- 
oped a DFT code 8 and applied it to the excitation and ionization energies of 
the atoms of the first transition series. The Becke exchange functional 9 cou- 
pled with the Lee- Yang-Parr correlation functional, 10 the (B-LYP) functional, 
gave exceptional results for the ionization energies (mean errors of ~ 0.1 eV) 
and acceptable results for excitation energies (mean errors of ~ 0.5eV). We 
were therefore interested in the predictions of gradient-corrected DFT for a 
"simple" series of ionic complexes, the transition metal fluorides: ScF 3 , TiF 4 , 
VF5, and CrFg. Experimental geometries and bond energies are in hand for 
the first three members of this series, thereby allowing a study of the accu- 
racy of the results as a function of basis set and functional, while the last 
one, CrF 6 , has been the subject of considerable controversy in recent years. 
On the experimental side, the debate centers about whether the vibrational 
spectrum observed by Hope, et al., 11,12 and used by them to assign Oh sym- 
metry to CrF 6 , indeed corresponds to CrF 6 , or is actually the signature of 
C1F5. 13,14 The debate among theorists centers around the relative stability 
of the Oh isomer relative to a trigonal prismatic geometry of D 3 h symme- 
try. 15-19 Early calculations preferred the D 3h structure but the most recent 
and extensive calculations suggest that the Oh isomer is the more stable. 

We present here the results of calculations on this series using both a lo- 
cal density approximation(LDA) and the gradient-corrected B-LYP approx- 
imation. The local density approximation comprises the Slater exchange 
functional and the Vosko-Wilks-Nusair parameterization of the electron gas 
correlation energy (S-VWN). Hartree-Fock calculations were also performed 
for comparision. Our results show that the B-LYP approximation yields rea- 
sonable binding energies for the transition metal fluorides using modest basis 
sets, but tends to overbind by 8-10 kcal/mol per bond. We have also inves- 
tigated a hybrid HF/DFT method due to Becke (the Becke3LYP method), 29 
and have found this to yield very good results (errors of the order of ~ 2 
kcal/mol per bond). The calculations support the claims that the octahedral 
structure of CrF 6 is of lower energy than the trigonal prismatic structure, 
our best calculation (Becke3LYP) favoring the octahedral structure by 14 
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kcal/mol. 



2 Method 

The DFT calculations were carried out using the MESA suite of programs, 20 
into which we recently added DFT modules. 8 The Kohn-Sham procedure 
used is strictly analogous to traditional Hartree-Fock calculations, with the 
replacement of the analytical exchange matrix with an exchange-correlation 
matrix calculated by numerical integration. 1 Analytic gradients of the energy 
were computed including the derivatives of the quadrature weights. 1 

The "Standard Grid-1" of Gill et al. 21 is used on the fluorines; on the met- 
als we use a simple grid of 100 concentric spheres, with radii chosen according 
to an Euler-MacLaurin scheme, 22 on which are positioned quadrature points 
chosen according to the scheme of Lebedev, 23 using the 194 point (23rd or- 
der) formulae. This 100-point radial grid for the transition metal has been 
shown to provide atomic energies accurate to 10~ 4 a.u. in a previous study. 8 
It is much larger than the fluorine grids and could be made significantly 
more economical with little loss in accuracy by applying a "grid pruning" 
technique. 21 ' 24 ' 25 

Most atomic calculations were performed using a spin-restricted open- 
shell procedure as described in a previous work. 8 

The hybrid HF/DFT, Becke3LYP, calculations were performed with the 
Gaussian 92/DFT programs. 30 A 100 point radial grid and 194 point angular 
grid was used on each atom. Bond energies were obtained using atomic 
energies calculated with a spin-unrestricted Kohn-Sham procedure. 

2.1 Basis sets 

Our initial studies on these molecules used a general contraction of the 
(lAs9p5d) Cartesian-Gaussian basis set of Wachters 26 augmented with the 
diffuse d function of Hay. 27 As described previously, 8 we have used Wachters' 
Hartree-Fock coefficients in our contraction, and contract the inner parts of 
the Is, 2s, 3s, 2p, 3p, and 3c? orbitals to obtain a [6s5p3rf] set using the gen- 
eral contraction scheme of Raffenetti. 28 The more diffuse primitives of each 
space are left free. This contracted set gave B-LYP total energies for Sc-Cr 
to within 10 -3 a.u. of the completely uncontracted results, and excitation 
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and ionization energies to within O.OleV. 8 With this metal basis set, we then 
examined various levels for the fluorine basis set in CrFg. Table 1 reports 
the optimal geometrical parameters, average bond energies, and relative sta- 
bility for the Oh and D 3 h conformers of CrF 6 in the B-LYP approximation. 
Examination of Table 1 reveals that inclusion of the d polarization func- 
tion (6-31G*) decreases the optimum bond length in both the Oh and D 3 h 
conformers by some 0.02 A. Consistent with this shorter bond length, the 
average bond energy increases by ~ 5 kcal/mol in both conformers. The 
relative energy is only slightly changed, the Oh conformer being more stable 
by ~ 9 kcal/mol. Because these species are expected to be ionic in character, 
calculations were performed in which diffuse s and p functions were added 
to the fluorine atom (6-31G+). Relative to the double-zeta 6-31G basis, this 
modification increased the bond lengths by 0.01 A, and decreased the average 
bond energy by ~ 5 kcal/mol. The addition of the diffuse functions had a 
more pronounced influence on the energy difference than did inclusion of the 
polarization function; the Oh form is 16.6 kcal/mol more stable than the D%h 
form in the (6-31G+) basis. A final calculation in which both diffuse and 
polarization functions were added to the fluorines (6-31G+*) yields a nearly 
identical relative stability of 16.3 kcal/mol. In the calculations to follow, 
then, the (6-31G+*) basis is used for the ligands. 

There were additional issues regarding the transition metal basis set which 
we wished to address before examining the series. First of all, folklore has it 
that Wachters' primitive set optimized for the atoms is sometimes inadequate 
in the valence 4s/p region in a molecule. Specifically, the primitive s space 
may be too diffuse, and the functions provided by Wachters' to describe the 4p 
atomic orbital are also more diffuse than those which might be needed to act 
as a polarization function in a molecule. Preliminary tests were therefore run 
on ScF 3 and TiF 4 in which the two most diffuse of Wachters' s primitives were 
replaced by tighter functions chosen to make the set more even-tempered. 
For Sc the new exponents were .12 and .04, for Ti they were .2 and .07. For 
ScF 3 the change in total energy with the modification of the s set was only 
~ 0.25 kcal/mol. Addition of a diffuse p function, a = 0.15, changed the 
total energy of ScF 3 by ~ 1 kcal/mol. Similar results were found for TiF 4 . 
Modification of the metal s and p set thus appears to be a minor effect in 
this series, presumably because of the ionic character of the bonding. 

A second issue concerns the relative importance of / functions in DFT cal- 
culations. Previous research using CI and coupled-cluster techniques discov- 
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ered a profound effect on the relative stability when / polarization functions 
were included on the metal. We found that addition of a metal / function, 
optimized in the B-LYP calculation, had a relatively minor effect on the en- 
ergy. In both ScF 3 and TiF 4 , addition of an optimized / function changed 
the total energy by ~ 5 kcal/mol, or ~ 1 kcal/mol in the average bond 
energy. In the discussion to follow, we compare results using the [6s5p3rf] 
contraction of Wachters' original basis and this basis augmented with an / 
function optimized in the B-LYP calculation [6s5p3dlf]. These optimized / 
exponents were found to be 0.4 for Sc, and 0.5 for the other three molecules. 

3 Results and Discussion 

3.1 ScF 3 , TiF 4 , and VF 5 

Table 2 compares theoretical and experimental bond lengths for ScF 3 , TiF 4 
and VF5. 31 The two major headings refer to the results using the [6s5p3rf] 
and [6s5p3dlf] metal atom bases. It can be seen that the additional metal 
/ function has little effect on the equilibrium distances, regardless of the 
theoretical method used. If any trend can be observed, it is a tendency for the 
/ function to decrease the bond lengths slightly (~ 0.01 A). Concentrating on 
the [6s5p3dlf] results, HF theory and the S-V WN aproximation are generally 
in good agreement with experiment. An exception occurs for ScF 3 where HF 
theory underestimates the bond length by 0.05A, with S-VWN being 0.08A 
too short. This exception was also noted in the S-VWN calculations of Sosa 
et al. 32 Although the basis sets, grids, and other details of our calculations 
differ from theirs, the agreement between the present S-VWN results and 
this earlier research is excellent. 

Perhaps most disappointing in Table 2 is the tendency of the B-LYP 
approximation to overestimate the bond lengths. For TiF 4 and VF 5 , where 
HF and S-VWN are in good agreement with experiment, the B-LYP results 
are ~ 0.04 — 0.05A longer than experiment. Similar behavior was obtained 
with B-LYP in a series of molecules comprised of first row atoms, 1 although 
the error was smaller, of the order of 0.02A. In the transition-metal fluoride 
series studied here, the origin of this behavior is associated with the Becke 
exchange functional. This can be seen from Table 3 where the equilibrium 
distances predicted by various functional combinations are presented. Note 
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that both the B-VWN and B-LYP approximations give distances which are 
too long, while S-VWN and S-LYP give distances in good agreement with 
experiment. The reason for this behavior is unknown; as is shown below, the 
hybrid HF/DFT version suggested by Becke largely remedies the problem. 29 

Table 4 compares average bond energies, E, for the various methods with 
experiment. The experimental values in Table 4 were obtained from en- 
thalpies of formation at 298.15°K. 33 The theoretical quantities reported are 
obtained by subtracting the energy of the molecule from the sum of the en- 
ergies of the atoms, and dividing by the number of fluorines. In order to 
compare the theoretical results directly with experiment, they should be cor- 
rected for the zero-point contribution and for the difference in enthalpies of 
formation between 298 and 0°K. In the case of TiF 4 , this thermodynamic in- 
formation is available from experiment. 34 To the experimental heat of atom- 
ization at 298°K, the zero-point contribution and differences in enthalpy for 
0°K, +5.5 kcal/mol and -3.5 kcal/mole, respectively, should be added to ob- 
tain a "corrected" experimental value. This "corrected" value thus increases 
by 2.0 kcal/mol, or 0.5 kcal/mol/bond for TiF4. As a result the average 
bond energy would be 140.4 instead of 139.9 kcal/mol, which can be com- 
pared directly with the theoretical value of 149.5 kcal/mol for BLYP. While 
vibrational frequencies are known for VF5, the other information needed to 
make corrections for the other molecules is not available. We estimate that 
these corrections would range from 0.3 to 0.7 kcal/mol for the average bond 
energies as one goes from ScF 3 to CrF 6 . These corrected estimates appear 
parenthetically in Table 4 

The first point to notice is that the /-functions have a minor influence 
on the average bond energy. They contribute ~ 2 kcal/mol per bond in the 
DFT calculations, and somewhat more for the HF calculation, up to ~ 4 
kcal/mol in VF 5 . As expected, the HF bond energies are in poor agree- 
ment with experiment. The error per bond grows quickly as the number of 
bonds increases; the error is ~ 60 kcal/mol per bond in VF 5 . This is shown 
graphically in Figure 1. On the other hand, the error in the S-VWN DFT 
approximation is more constant as the size of the molecule increases. The 
S-VWN approximation overbinds the molecule, by 30, 33, and 39 kcal/mol 
per bond in ScF 3 , TiF 4 , and VF 5 , respectively. The B-LYP bond energies are 
in much better agreement with experiment. Like S-VWN, this approxima- 
tion also overestimates the bond energy, but by only 8, 9, and 13 kcal/mol, 
respectively. This error is somewhat larger than the typical error observed in 



7 



studies of organics, but it should be kept in mind that quite extensive theoret- 
ical treatments including very large one- and many-electron bases would be 
necessary in order to obtain this degree of accuracy in conventional HF-based 
correlation schemes. 

It is sometimes argued that more accurate bond energies are obtained 
when using spin-unrestricted Kohn-Sham (UKS) atomic results, as opposed 
to the spin-restricted open-shell method used here. Since the DFT approx- 
imations generally overbind the molecule, and the UKS method will always 
yield atomic energies identical to or lower than the ROKS method, the use 
of UKS atomic energies will certainly decrease the bond energies, thereby 
improving the apparent agreement with experiment. In order to investigate 
the magnitude of this effect, we calculated UKS atomic energies using Gaus- 
sian92/DFT. 30 In the [6s5p3d]/[6 - 31G + *] basis set, the average bond 
energies were 153.4, 146.8 and 122.3 for ScF 3 , TiF 4 and VF 5 . The difference 
between our ROKS and these UKS results is small (~ 1 kcal/mole per bond). 

Becke 29 has argued that the residual tendency of the gradient-corrected 
DFT methods to overbind stems from the exchange functional and has intro- 
duced a hybrid method in which a component of the "exact" Hartree-Fock 
exchange is retained in the energy expression. We have therefore also exam- 
ined a variant of this hybrid approach, the Becke3LYP method. 

The [6s5p3dl/]/[6 — 31G + *] basis set was used, giving the results shown 
in Table 5. The results are in very good agreement with experiment: the 
difference in the average bond energies between the raw Becke3LYP results 
and the "corrected" experimental values is 2.0, 1.4 and 1.7 kcal/mol per 
bond. Bond lengths are also significantly improved over the B-LYP approach. 
They are approximately 0.01A too long, except in the case of SCF3, where the 
bond length is underestimated by 0.05A. Given the success of the method 
for the other molecules, this experimental determination should probably be 
revisited. 

3.2 CrF 6 

Before discussing the present results for CrF 6 , a brief review of the history of 
the problem is useful. In 1963, Hellberg, Muller and Glemser 36 reported the 
preparation of a volatile yellow substance which they identified as CrF 6 . On 
the basis of vibrational data, Hope, et al. 11 assigned octahedral (Oh) sym- 
metry to the molecule. The first ab initio theoretical treatment of d° MLq 
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complexes was reported by Kang, Albright and Eisenstein in 1989 15 ; these 
studies predicted octahedral geometry for CrFg at the MP2/minimum basis 
set level of theory. This prediction was countered by the calculations of Mars- 
den and Wolynec 16 who concluded that CrF 6 was probably trigonal prismatic 
{D 3 h) based on the predictions of coupled-cluster theory with double substi- 
tutions (CCD) in a double-zeta basis set. Jacob and Willner 13 carried out 
new experiments and argued that CrF 6 had not yet been prepared, and that 
earlier experiments had merely produced CrF 5 . Hope, Levason and Ogden 12 
responded with a more detailed analysis of their data and again concluded 
that CrFg had octahedral geometry. 

More extensive theoretical treatments then came into play. Pierloot and 
Roos 17 performed CASPT2 (complete active space with second-order M0ller- 
Plesset perturbation theory) calculations in a large basis set, and found that 
the octahedral isomer was 50 kcal/mol lower in energy than the prismatic 
isomer. Neuhaus, Frenking, Huber and Gauss 18 discovered two interesting 
facts. The first is a large differential effect in the SCF energies when / func- 
tions are included in the Cr basis set. SCF calculations without / functions 
showed the prismatic isomer to be preferred, but with / functions the oc- 
tahedral isomer was favored by 11 kcal/mol. Neuhas, et al. then explored 
a hierarchy of coupled cluster expansions with a smaller double-zeta basis 
and discovered a second interesting point: a large differential effect in the 
correlation energies of the two isomers associated with the contribution of 
triple excitations. In their double-zeta basis it was nearly as large as the 
differential effect of / functions at the SCF level. By combining these two 
results, Neuhas et al. concluded that the octahedral isomer was lower in en- 
ergy by some 20 kcal/mol. Meanwhile, back in the lab, a 1992 experimental 
paper from Jacobs, Miiller, Willner, Jacob and Burger 14 once again argued 
that CrF 6 had not in fact been prepared. Most recently, in 1994, Marsden, 
Moncrieff and Quelch 19 reported correlated calculations using CISD, CCSD 
and CCSD(T) methods with a large basis including / functions on the metal. 
Their work confirmed the sensitivity of the results to metal / functions and 
triple excitations. Their most extensive calculation favored the octahedral 
isomer by 14.2 kcal/mol. Thus, at the present time, theory is converging to 
the position that the octahedral form is more stable than the prismatic. The 
magnitude of the difference is still not settled, particularly given the large 
effect of triple excitations on the relative stability. 

The results of the B-LYP calculations for CrFg are presented in Table 2 
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and Table 4, and those using the Becke3LYP approximation in Table 5. 
Once again, the Becke3LYP method yields shorter bond lengths and smaller 
binding energies than B-LYP. The Becke3LYP equilibrium bond distances for 
both the Oh and forms are similar, 1.73A and 1.74 A respectively. The 
optimum bond angle in the D 3 h isomer, i.e. that between a three-fold axis 
and a Cr-F bond, was found to be 50.45°, in good agreement with earlier 
HF calculations which found 50.5°. The average bond energy follows the 
trend established in the earlier members of the series of decreasing strength 
with increasing number of bonds. The difference in energy between the two 
isomers is 13.8 kcal/mol. This result is very close to the latest CCSD(T) 
calculations of Marsden, et al. who found 14.2 kcal/mol. 

Since the experimental efforts at identifying the molecule have focussed 
on the vibrational spectrum, we present in Table 6 the B-LYP harmonic 
frequencies of the Oh isomer. The two infrared active modes occur at 329 
and 717cm" 1 . The former is in excellent agreement with the experimental 
bend observed 12 at 332cm -1 . The latter is in reasonable agreement with 
the intense mode observed at 760cm" 1 . Although there are not many cases 
reported yet, it appears that the gradient-corrected approximations tend to 
underestimate metal-ligand stretching frequencies. This would be consistent 
with the observed tendency of B-LYP to overestimate bond lengths. As dis- 
cussed above, the Becke3LYP approximation significantly improved both the 
bond lengths and the bond energies relative to B-LYP; we therefore also ex- 
amined the vibrational frequencies predicted by Becke3LYP (Table 6). The 
theoretical predictions of 768 and 336 cm -1 are in excellent agreement with 
the experimentally observed bands. Furthermore, the ratio of the IR intensi- 
ties of the Becke3LYP modes is 0.081, which agrees with the corresponding 
ratio in the experiments of Jacob and Willner, 13 0.080. These calculations 
strongly support both the conclusion that CrF 6 has been prepared and the 
assignment of Oh symmetry. 

4 Conclusions 

As has been noted many times before, 37 an advantage of the DFT approach 
is that the density, and hence the energy, converges much more rapidly with 
basis set than does the correlation energy in configuration-interaction or 
coupled-cluster approaches. A case in point in the present work is the rela- 
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tive insensitivity of the results to the presence of / functions on the metal. 
This can be contrasted with the research reviewed above. In Hartree-Fock 
based correlation approaches the / functions must provide for angular corre- 
lation among the d electrons, while in DFT calculations they play the role of 
polarization functions only. The small effect they have on the structure and 
bond energies in this series is presumably a reflection of the ionic character 
of the bonding. In the ionic limit, the metal density is roughly spherical, and 
the dominant polarization effects should come from p functions. 

The bond energies provided by the B-LYP approximation are in fairly 
good agreement with experiment. The decrease in average bond energy 
down the series is faithfully reproduced, and the absolute values are within 
~ 10 kcal/mol of experiment even in this relatively modest double-zeta- 
polarization basis set. The B-LYP distances, however, appear to be consis- 
tently too long by ~ 0.04 A, with the exception of ScF 3 . The anomalous 
behavior of ScF 3 suggests a problem with the experimental data. 

Both the bond distances and energies are greatly improved in the hybrid 
HF/DFT Becke3LYP approach. Distances are within O.OlA and average 
bond energies within 2 kcal/mole per bond of experiment. This success is 
impressive, especially considering that the molecules studied here are quite 
different from those used by Becke to generate his fit. We note that Ricca 
and Bauschlicher 7 have recently found that the Becke3LYP approximation 
provides excellent thermochemistries in the series Fe(CO)^. The bonding in 
these carbonyls is certainly quite different from the ionic limit studied in the 
current contribution; the approach appears very promising given the small 
sample of transition metal complexes studied thus far. 

Finally, our B-LYP calculations suggest that the octahedral form of CrFg 
is preferred over the D%h form by ~ 17 kcal/mol; this number is reduced 
to 14 kcal/mol in the Becke3LYP approximation. The Becke3LYP IR-active 
vibrational frequencies and relative intensities computed for the Oh struc- 
ture are in excellent agreement with the experimental infrared spectrum and 
support the original assignment of Hope, et al. 
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Basis 


6-31G 


6-31G* 


6-31G+ 


6-31G+* 


r(O h ) 


1.78 


1.76 


1.79 


1.75 


r(D 3h ) 


1.79 


1.77 


1.80 


1.78 


E{O h ) 


93.1 


98.3 


88.2 


92.6 


E(D 3h ) 


91.5 


96.8 


85.4 


89.9 


E(D 3h -O h ) 


9.5 


9.0 


16.6 


16.3 



Table 1: Equilibrium bond lengths, average bond energies, and conformer 
stability for CrF 6 . The columns denote various F bases used in conjunction 
with the [6s5p3<i] contraction of the Wachters-Hay primitive set for Cr (r in 
A, E in kcal/mol). 
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[6s5p3d 






[6s5p3dlf] 


Expt 


Molecule 


HF 


S-VWN 


B-LYP 


HF 


S-VWN 


B-LYP 




ScF 3 


1.860 


1.831 


1.868 


1.858 


1.829 


1.865 


1.91 


TiF 4 


1.747 


1.745 


1.781 


1.743 


1.740 


1.776 


1.745 


VF 5 r ax 


1.728 


1.747 


1.786 


1.718 


1.738 


1.778 


1.734 




1.685 


1.709 


1.748 


1.681 


1.704 


1.744 


1.703 


CrF 6 O h 






1.753 






1.765 




D 3h 






1.782 






1.775 





Table 2: Optimum bond lengths(A) from S-VWN, B-LYP and HF calcula- 
tions on the transition metal fluoride series. The calculations were carried 
out with the 6-31G+* basis on the fluorine atoms and either the [6s5p3rf] 
or the [6s5p3dlf] contraction described in the text. The experimental re- 
sults are summarized in Ref. 31. For CrF 6 (D 3 / l ), the optimum bond angle 
between a three-fold rotation axis and the Cr-F bond is 50.44° and 50.45° in 
the [6s5p3ci] and [Qs5p3dlf] metal bases, respectively. 



[6s5p3d] Expt 



Molecule 


S-VWN 


S-LYP 


B-VWN 


B-LYP 




ScF 3 


1.831 


1.826 


1.874 


1.868 


1.91 


TiF 4 


1.745 


1.739 


1.787 


1.781 


1.745 


VF 5 r ax 


1.728 


1.740 


1.794 


1.786 


1.734 


f'eq 


1.685 


1.703 


1.756 


1.748 


1.703 



Table 3: Optimum bond lengths(A) from S-VWN, S-LYP, B-VWN, and B- 
LYP calculations on the transition metal fluoride series. The calculations 
were carried out with the 6-31G+* basis on the fluorine atoms and the 
[6s5p3rf] contraction described in the text. The experimental results are 
summarized in Ref. 31. 
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[6s5p3d] [6s5p3dlf] Expt 



Molecule 


HF 


S-VWN 


B-LYP 


HF 


S-VWN 


B-LYP 




ScF 3 


133.8 


175.6 


154.4 


135.7 


177.3 


155.6 


147.4(147.7) 


TiF 4 


91.2 


171.5 


147.8 


94.1 


173.6 


149.5 


139.9(140.4) 


VF 5 


49.4 


148.8 


123.6 


53.3 


151.2 


125.5 


112.1(112.7) 


CrF 6 O h 






92.6 






94.5 




D 3h 






89.9 






91.7 





Table 4: Average bond energies £(kcal/mol) from HF, S-VWN, and B-LYP 
calculations on the transition metal fluoride series. The calculations were 
carried out with the 6-31G+* basis on the fluorine atoms and either the 
[6s5p3rf] or [6s5p3dlf] contraction described in the text. The experimental 
values are obtained from heats of formation at 298°K. 33 Estimates of the 
experimental values "corrected" for zero-point energy and enthalpy differ- 
ences between and 298°K are given parenthetically (see text); these may 
be compared directly to the theoretical results. 



Becke3LYP Experiment 



Molecule 




E 


'"eg 


E 


ScF 3 




1.852 


149.7 


1.91 


147.4(147.7) 


TiF 4 




1.756 


141.8 


1.745 


139.8(140.4) 


VF 5 


fax 


1.750 


114.4 


1.734 


112.1(112.7) 




T'eq 


1.717 




1.703 




CrF 6 


o h 


1.732 


78.8 






CrF 6 


D 3h 


1.743 


76.5 







Table 5: Equilibrium bond lengths and average bond energies calculated with 
Becke3LYP method compared to experimental values. The [6s5p3dlf}/6- 
31G+* basis set was used. Bond energies were calculated using UKS energies 
for atoms. The experimental values are obtained from Ref. 31 and Ref. 33; 
see caption for Table 4 and text. 
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B-LYP 


Becke3LYP 


Experiment 


Mode 


We(cm _1 ) 


w e (cm _1 ) 




hu 


157.8 


145.2 




tlu 


329.0 


335.7(27.06) 


332(0.08) 


hg 


344.2 


363.3 




e 9 


554.7 


591.8 




dig 


648.8 


708.4 




tlu 


717.0 


767.7(333.67) 


760(1.0) 



Table 6: Vibrational frequencies in CrF 6 . The calculations were carried out 
with the 6-31G+* basis on the fluorine atoms and the [6s5p3dlf] metal 
contraction described in the text. They were obtained by finite differencing 
analytic gradients at the Oh equilibrium geometry. Numbers in parentheses 
are the IR intensities in KM/mol. The experimental frequencies and relative 
intensities are from Ref. 13. 
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6 Figure Captions 
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□ 














o 










o 
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Figure 1: The average bond energy in kcal/mol for the transition metal 
fluoride series for HF(O), S-VWN (+), and B-LYP (□), Becke3LYP (x) and 
experiment (A). 
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